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ABSTRACT 

Broad Absorption Line (BAL) QSOs have been suggested to be youthful 
super-accretors based on their powerful radiatively driven absorbing outflows and 
often reddened continua. To test this hypothesis, we observed near IR spectra 
of the H/3 region for 11 bright BAL QSOs at redshift z ~2. We measured these 
and literature spectra for 6 BAL QSOs, 13 radio-loud and 7 radio-quiet non-BAL 
QSOs. Using the luminosity and H/3 broad line width to derive black hole mass 
and accretion rate, we find that both BAL and non-BAL QSOs at z ~2 tend 
to have higher L/L-^dd than those at low z - probably a result of selecting the 
brightest QSOs. However, we find that the high z QSOs, in particular the BAL 
QSOs, have extremely strong Fe II and very weak [O III], extending the inverse 
relationship found for low z QSOs. This suggests that, even while radiating near 
LEdd, the BAL QSOs have a more plentiful fuel supply than non-BAL QSOs. 
Comparison with low z QSOs shows for the first time that the inverse Fe II - 
[O III] relationship is indeed related to L/L^d, rather than black hole mass. 

Subject headings: quasars: emission lines — quasars: general — galaxies: active 
- galaxies: evolution — accretion, accretion disks — black hole physics 



1. INTRODUCTION 

Broad Absorption Line (BAL) QSOs show very broad blueshifted absorption troughs 
indicating outflows up to ~ 0.1c, with covering fractions > 0.2 (Hewett & Foltz 2003), and 
soft X-ray hydrogen column densities often exceeding 10 23 cm -2 (Gallagher et al. 2002). 
These outflows carry a significant fraction of the accretion power and there is evidence that 
they are radiatively driven (Arav 1997), suggesting that BAL QSO black holes accrete close 
to the Eddington limit. Many BAL QSOs, especially those with low ionization absorption 
troughs, are reddened (Low et al. 1988; Becker et al. 2000; Hall et al. 2002), suggesting 
that BAL QSOs may be young QSOs ejecting a birth-cocoon of dusty gas to reveal the 
unobscured optical QSO (Fabian 1999; Haehnelt, Natarajan & Rees 1998). With the 
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discovery of many BAL QSOs at z ~ 2, we are able to explore this idea by comparing 
Eddington accretion ratio (L/L Edd ) directly with fueling indicators in BAL and non-BAL 
QSOs. 

We can estimate M B h and hence L/L E dd for z ~ 2 QSOs using observed continuum 
luminosity and FWHM Wf3 in the following way: If QSO Broad Line Region (BLR) gas 
motion is dominated by gravity, the black hole mass (Mbh) can be calculated from the 
typical orbital radius of BLR clouds and their typical velocity: Mbh — v 2 r/G. BLR radius r 
can be represented by BLR size (-Rblr) derived from Wf3 reverberation mapping: -Rblr can 
be estimated from continuum luminosity at rest wavelength 5100A: Rblr oc La(5100A) 066 
(Peterson & Wandel 1999; Kaspi et al. 2000; Vestergaard 2002); v can be represented by 
H/3 full width at half maximum (FWHM). The M B h so derived agree with the M B h derived 
independently from the Mbh - bulge luminosity and Mbh oc <t 4 relationships for nearby 
QSOs, where a is the velocity dispersion of the bulge (Laor 1998; Gebhardt et al. 2000; 
McLure & Dunlop 2002). 

The fueling could be indicated by correlated emission line properties. One of the 
strongest sets of relationships among QSO properties describes the increasing strength of 
optical Fe II with decreasing [O III] A5007, increasing steepness of the soft X-ray continuum, 
and decreasing width and stronger blue wing of the broad H/3 emission line (Grupe et al. 
1999; Shang et al. 2003; Boroson & Green 1992, hereinafter BG92). Here we call this 
set of relationships Boroson & Green Eigenvector 1 (BGEV1). It has been suggested that 
BGEV1 is driven by L/L Edd for the following reasons. 

• For a given luminosity, narrower Wf3 corresponds to smaller Mbh and higher L/L Edd . 

• By analogy with Galactic black hole binary systems, it has been suggested that steeper 
X-ray spectra indicate higher L/L Edd (Pounds, Done & Osborne 1995; Done et al. 
1995). Near-Eddington accretion results in geometrically thick accretion disks which 
produce excess soft X-ray photons (e.g. Tanaka & Shibazaki 1996). 

• Strong optical Fe II emission and weak [O III] emission indicate high optical depth, 
high density nuclear gas that may provide an abundant fuel supply for near-Eddington 
accretion. Optically thick gas reduces ionizing photons that reach the [O III] narrow 
line region (e.g. Boroson & Green 1992). Strong Fe II suggests that this gas may 
have been metal-enriched by star formation. The few known BAL QSOs at low z tend 
to have narrow H/5, weak [O III] and strong Fe II emission (Turnshek et al. 1997; 
Weymann et al. 1991) suggesting high L/L Edd and an abundant fuel supply. 

Most BALs are discovered in QSOs with z > 1.5, where the broad C IV troughs are 
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shifted into the optical window. For those redshifts the B./3 region is shifted to the near 
infrared. We therefore performed near infrared spectroscopy, taking advantage of new QSO 
surveys that do not use UV selection and yield many more BAL QSOs than previously 
recognized (Becker et al. 2000; Reichard et al. 2003). In this letter, we investigate whether 
BAL QSOs have extreme BGEV1 properties, using FWHM H(3 and continuum luminosity 
to calculate L/Ls&d, and then compare L/L^dd with Fe II and [O III] strengths. Detailed 
discussion of the observations, reductions and analysis, with results from a larger data set, 
will be published later. 



2. SAMPLES AND DATA 

2.1. Samples and Observations 

Our high redshift (z ~ 2) sample consists of 37 QSOs. We observed 11 BAL QSOs 
from the Large Bright QSO Survey, the FIRST Bright QSO Survey and Weymann et al. 
(1991). QSOs of similar redshift and luminosity were observed by Mcintosh et al. (1999a, 
hereinafter M99a). We included their observations having acceptable signal-to- noise ratio, to 
make a total sample of 17 BAL QSOs, 13 radio-loud (RL) QSOs and 7 non-BAL radio-quiet 
(RQ) QSOs. We and M99a chose the brightest optically-selected QSOs at z ~ 2. 

Our new spectra were obtained with the CGS4 near infrared spectrograph on the United 
Kingdom Infra- Red Telescope (UKIRT). The detector was the 256x256 InSb NICMOS array. 
With the 40 1 mm -1 low resolution grating, a slit width of 1.2" yielded a resolution of 3 
pixels as measured from Ar and Xe comparison lamp spectra (~ 600 km s _1 in J band and 
~ 750 km s _1 in H band). 

The two-dimensional spectral images were reduced using UKIRT's ORAC data reduction 
pipeline. We used standard IRAF tasks to optimally extract the spectra and calibrate their 
wavelength scales. Atmospheric features were removed and spectral shape calibration was 
done by division by the observed spectra of F-G stars using a temperature implied by the 
spectral type. The F-G stars and the wavelength comparison lamps were observed within 
0.1 airmasses of the QSOs. Our reduced spectra have signal-to-noise ratios of from 10 to 20 
per 3-pixel bin near LL3. 
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2.2. Spectral Decomposition 

Measuring line and continuum properties in the B./3 region is complicated by the blend- 
ing of many broad optical Fe II lines. We used SPECFIT (Kriss 1994) within IRAF to 
deblend spectral components for both our UKIRT and the M99a spectra. We decomposed 
each spectrum into a powerlaw continuum, broad and narrow H/3 components, a broad H7 
component, a narrow-line [O III] doublet and Fe II emission blends. 

The fixed parameters were the power-law continuum parameters (slope and normaliza- 
tion estimated by eye), the [O III] AA4959,5007 doublet ratio of 2.94, the H7/H/3 ratio of 
0.36, and the narrow H/3 to [O III] A5007 ratio of 0.1 (Veilleux & Osterbrock 1987). Except 
for some M99a spectra with known broad [O III] lines, the narrow H/3 and [O III] lines were 
represented by single Gaussians with width equal to the instrumental resolution (BG92). 
The Fe II emission blends were represented by the BG92 I Zw 1 template. All broad lines 
including Fe II were assumed to have the same Gaussian profile. Rest wavelength ratios were 
constrained except for some Mcintosh et al. objects with shifted [O III] lines (Mcintosh et 
al. 1999b). Free parameters were the intensities of [O III], broad H/3, the Fe II blends, and 
broad line widths. 

In many cases, since a single Gaussian profile was inadequate to represent the broad 
H/9 line, we adopted the following procedure after running SPECFIT. We first subtracted 
all fitted components except broad H/3 and smoothed the remaining H[3 by fitting multiple 
Gaussian profiles, each at least as wide as the instrumental resolution, then measured the 
integrated flux and FWHM of the smoothed profile. We tested our measurement method 
on BG92 spectra, finding no systematic difference from their results. For objects without 
measurable [O III] or Fe II lines, we estimated 3a upper limits based on the signal-to- 
noise ratio and spectral resolution of each spectrum. The typical rms uncertainties of our 
measurements for the high redshift QSOs are 10% for [O III] A5007 equivalent width (EW) 
for objects with detected [O III], 10% for FWHM H/3, and about 20% for Fe II/H/3 ratios. 
Figure 1 shows the fits for three typical spectra. The BAL QSO spectrum is from our new 
UKIRT data. The RL and RQ QSO spectra are from our refitting of the M99a data. 
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Fig. 1. — Typical observed spectra (solid lines) and fitted models (dotted lines) for high 
redshift QSOs. The underlying solid line is the fitted continuum. 
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2.3. Continuum Luminosities 

To calculate Mbh and L/L^di we need L\ at rest frame wavelengths of 5100A and 
3000A to estimate the BLR size, and the bolometric luminosity: L^ Q \ = 8.3 • ALa(3000A) 
(Laor 2000). BAL QSOs tend to have reddened continua compared with optically selected 
non-BAL QSOs, so the L x were calculated using H band magnitudes from the the 2MASS 
survey 1 , assuming f v oc z/ -0 5 . Using the most recent BAL QSO colors derived from the Sloan 
Digital Sky Survey (Tolea, Krolik & Tsvetanov 2002), we estimate a generous upper limit 
for extinction at the observed H band correponding to ~0.8 mag for z ~ 2, corresponding to 
under-estimates of black hole mass and Eddington ratio by factors of < 10~°- 2 and < 10~ ai 
respectively. Such small corrections would not change our results significantly. We assumed 
a cosmological model with H Q = 70 km s _1 Mpc" 1 , Q m = 0.3 and = 0.7 (Freedman 
2002). Table 1 lists the line measurements and calculated L/L^a for our high z QSOs. The 
complete table is available electronically. 



htt p : / / www. ipac .caltech . edu / 2mass / 
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Table 1. z ~ 2 QSO properties 



Name 


Type 


EW [0 III] 

A(rest) 


Fe II/R/3 


M BH 
1O 9 M 


L/ L Edd 


0046+0104 


BAL 


11.3 


3.0 


3.3 


0.78 


0052+0101 


RQ 


6.4 


5.6 


5.1 


0.9 


0052+0140 


RQ 


5.7 


2.9 


0.7 


3.2 


0126+2559 


RL 


20.9 


1.7 


0.4 


6.76 


0157+7442 


RL 


19.7 


1.3 


19.4 


0.42 


0228-0337 


RL 


23.9 


1.2 


1.9 


1.07 


0228-1011 


BAL 


7.0 


3.3 


9.5 


0.26 


0424+0204 


RL 


50.6 


1.3 


15.5 


0.30 


0427-1302 


RL 


29.9 


1.0 


8.6 


0.47 


0555+3948 


RL 


14.9 


<1.8 


1.9 


1.71 


0724+4159 


BAL 


<5.6 


7.4 


1.6 


1.91 


0841+7053 


RL 


4.8 


5.1 


4.0 


1.86 


0845+3420 


BAL 


5.6 


4.2 


16.4 


0.18 


0913+3944 


BAL 


<3.0 


8.7 


0.1 


9.48 


0934+3153 


BAL 


<1.8 


3.1 


5.1 


0.80 


1013+0851 


BAL 


<1.3 


4.6 


14.3 


0.39 


1054+2536 


BAL 


<2.7 


4.2 


2.5 


1.07 


1106-1821 


RQ 


12.9 


2.3 


4.4 


1.0 


1225+2235 


RQ 


8.9 


1.9 


35.8 


0.3 


1228+3128 


RL 


<0.9 


5.1 


40.8 


0.39 


1231+0725 


RL 


8.3 


<0.7 


3.5 


1.04 


1233+1304 


BAL 


<6.2 


3.1 


8.3 


0.41 


1234+1308 


BAL 


<3.0 


2.4 


2.8 


1.04 


1249-0559 


BAL 


<2.3 


5.0 


22.9 


0.51 


1250+2631 


RQ 


10.6 


1.5 


11.4 


1.2 


1311-0552 


BAL 


<1.3 


6.0 


4.8 


0.94 


1333+1649 


RL 


18.7 


1.7 


17.5 


0.58 


1348-0353 


RQ 


4.3 


3.3 


3.9 


1.3 


1418+0852 


RQ 


8.2 


2.3 


4.6 


1.0 


1420+2534 


BAL 


<4.9 


11.1 


3.0 


1.07 
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Table 1 — Continued 



Name 


Type 


EW [0 III] 


Fe II/H/3 


M BH 


L/ -^Edd 






A(rest) 




1O 9 M 




1436+6336 


RL 


11.6 


1.1 


16.4 


0.31 


1445+0129 


BAL 


<3.7 


2.6 


1.1 


1.13 


1445-0023 


BAL 


<4.2 


9.5 


8.0 


0.17 


1451-2329 


RL 


22.4 


1.4 


4.8 


1.20 


1516+0029 


BAL 


<2.1 


1.4 


6.7 


0.48 


2215-1744 


BAL 


20.3 


1.7 


15.4 


0.25 


2312+3847 


RL 


28.4 


5.3 


5.4 


0.63 
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3. RESULTS AND DISCUSSION 

3.1. Are BAL QSOs Extreme Accretors? 

Figure 2 plots FWHM Kf3 against bolometric luminosity. Lines of constant M B h and 
L/L Edd are indicated. High redshift QSOs (large symbols) are compared with the 87 low 
redshift BG92 QSOs from which the original BGEV1 relationships were derived (small sym- 
bols). Different sub-samples of QSOs are distinguished. The L/L Edd for our flat-spectrum 
RL QSOs might be overestimated if the continuum is beamed (Padovani & Urry 1992), 
and FWHM H.[3 may underestimate the true virial speed as a result of viewing disk motions 
pole-on (e.g. Krolik 2001). The brightest QSOs at z ~ 2 are the most luminous, introducing 
a bias toward the highest L/L Edd and M B h- The high z BAL QSOs, while having high 
L/L Edd , are not clearly different from the high z non-BAL QSOs, even after applying the 
orientation corrections. 
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Fig. 2. — FWHM H/9 versus bolometric luminosity. Lines of constant Mbh and constant 
L/L^dd are shown. See Fig. 3 for symbols. The arrows show possible corrections to the high 
z RL QSO data for velocity projection and continuum beaming. 



3.2. The Fe II versus [O III] Anti-correlation 



Figure 3 shows that high z QSOs, in particular the BAL QSOs, extend the inverse 
relationship between Fe II/H/3 and EW [O III] to stronger Fe II and weaker [O III]. Thirteen 
out of seventeen high z BAL QSOs do not even have detectable [O III] emission lines. Taking 
into account the [O III] upper limits, the generalized Kendall's tau correlation coefficient 
indicates that the 1-tailed probability for the Fe II versus [O III] correlation to arise by 
chance is < 0.01% for either the entire sample or the high z QSOs alone. The trend for low 
z RL and BAL QSOs to lie at opposite extremes has been noted before (BG92). The same 
appears to be true among the high z QSOs. The 1-tailed probability of high z BAL QSOs 
and low z RQ QSOs having the same line strength distribution is < 0.1% for Fe II/H/9 and 
< 0.1% for EW [O III]. 
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Fig. 3.— Intensity ratio Fe Il/H/3 versus EW [O III] A5007. 



3.3. Eddington Ratio Indicator 



. Figures 4 and 5 show Mbh and L/L Edd versus Fe II/H/3. In Fig. 4, the z ~ 2 
QSOs lie at significantly larger Mbh; in Fig. 5, these high redshift QSOs extend the low z 
relationship to higher L/L Edd . For both high and low z QSOs, the 2-tailed probability of the 
L/L Edd versus Fe II/H/3 correlation arising by chance is < 0.01%. The M B h versus Fe II/H/3 
correlation for low z QSOs may simply reflect the inverse dependence of Mbh on L/L Edd , 
given a limited luminosity range. There is also an inverse correlation between L/L Edd and 
EW [O III] for both high and low z QSOs (not shown), with most BAL QSOs lying at 
the weak [O III] extreme (see Fig. 3). The 2-tailed probability of this correlation arising 
from unrelated variables is < 0.01%. The above correlations demonstrate that BGEV1 is 
indeed related to L/L Edd rather than Mbh- If BGEV1 properties are L/L Edd indicators, why 
do high z BAL and non-BAL QSOs have different BGEV1 properties despite their similar 
L/L Edd ? One possibility is that all our high z QSOs belong to the same parent population 
but BAL QSOs are viewed from special angles. This would imply that BGEV1 emission line 
properties are affected by orientation. More likely, the BGEV1 emission line relationships 
may be only indirectly related to Lj ' L Edd . They may depend more directly on the availability 
of fuel around the black hole. Stronger Fe II and weaker [O III] correspond to an abundance 
of cold gas, which could fuel high Eddington ratio accretion. For low L/L Edd objects, which 
are found mostly in the low z sample, increasing the fuel supply increases the accretion 
rate. At higher luminosities (high z) BAL QSOs may be the youngest QSOs with the most 
abundant fuel supplies, as indicated by very strong Fe II and very weak [O III] emission, but 
there is not a corresponding increase in accretion ratio because they are unable to radiate 
at L > L Edd 
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Fig. 4. — M B h versus Fe II/H/3. See Fig. 3 for symbols. 
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Fig. 5. — L/L E dd versus Fe II/H/5. See Fig. 3 for symbols. 



-16- 



We thank P. Hirst for observing support, D. Wills for comments on a draft, D. H. 
Mcintosh for generously providing H band spectra for high redshift QSOs, T. A. Boroson for 
spectra of the BG92 sample, and Roc Cutri for providing unpublished 2MASS magnitudes. 
The United Kingdom Infra-Red Telescope is operated by the Joint Astronomy Centre on 
behalf of the U.K. Particle Physics and Astronomy Research Council. 

REFERENCES 

Arav, N. 1997, in ASP Conf. Ser. 128, Mass Ejection from Active Galactic Nuclei, ed. N. 
Aravm, I. Shlosman, & R. J. Weymann, 264 

Becker, R. H., White, R. L., Gregg, M. D., Brotherton, M. S., Laurent-Muehleisen, S. A. & 
Arav, N. 2000, ApJ, 538, 72 

Boroson, T. A. 2002, ApJ, 565, 78 

Boroson, T. A., & Green, R. F. 1992, ApJS, 80, 109 (BG92) 
Brandt, W. N., Laor, A. & Wills, B. J. 2000, ApJ, 528, 637 
Corbin, M. 1997, ApJS, 113, 245 

Done, C, Pounds, K. A., Nandra, K. & Fabian, A. C. 1995, MNRAS, 275, 417 
Fabian, A. C. 1999, MNRAS, 308, L39 
Freedman, W. 2002, Int. J. Mod. Phys., A17S1, 58 

Gallagher, S. C, Brandt, W. N., Chartas, G. & Garmire, G. P. 2002, ApJ, 567, 37 

Gebhardt, K., Kormendy, J., Ho, L. C. & Bender, R. 2000, ApJ, 543, 5 

Grupe, D., Beuermann, K., Mannheim, K. and Thomas, H.-C. 1999, A&A, 350, 805 

Haehnelt, M. G., Natarajan, P. & Rees, M. J. 1998, MNRAS, 300, 817 

Hall, P. B. et al. 2002, ApJS, 141, 267 

Hewett, P. C, & Foltz, C. B. 2003, AJ, 125, 1784 

Kaspi, S., Smith, P. S., Netzer, H., Maoz, D. Jannuzi, B. T. & Giveon, U. 2000, ApJ, 533, 
631 



-17- 

Kriss, G. 1994, in ASP Conf. Ser. 61, Astronomical Data Analysis Software and Systems III, 
ed. D. R. Crabtree, R.J. Hanisch, & J. Barnes, 437 

Krolik, J. H. 2001, ApJ, 551, 72 

Lacy, M., Laurent-Muehleisen, S. A., Ridgway, S. E., Becker, R. H. & White, R. L. ApJ, 
505, L83 

Laor, A., Fiore, F., Elvis, M., Wilkes, B. J. & McDowell, J. C. 1997, ApJ, 477, L93 
Laor, A. 1998, ApJ, 505, L83 
Laor, A. 2000, ApJ, 543, Llll 

Low, F. J., Cutri, R. M., Huchra, J. P. & Kleinmann, S. G. 1988, ApJ, 327, L41 

Mcintosh, D. H., Rieke, M. J., Rix, H.-W., Foltz, C. B., Weymann, R. J. 1999, ApJ, 514, 40 
(M99a) 

Mcintosh, D. H., Rix, H.-W., Rieke, M. J., Foltz, C. B. 1999, ApJ, 517, L73 (M99b) 

McLure, R. J. & Dunlop, J. S. 2002, MNRAS, 331, 795 

Padovani, P. & Urry, C. M. 1992, ApJ, 387, 449 

Peterson, B. M. & Wandel, A. 1999, ApJ, 521, L95 

Pounds, K. A., Done, C. & Osborne, J. P. 1995, MNRAS, 277, L5 

Reichard, T. A. et al. 2003, AJ, 125, 1711 

Shang, Z., Wills, B. J., Edward, L. R., Wills, D., Laor, A., Xie, B. & Yuan, J. 2003, ApJ, 
586, 52 

Storey, P. J. & Zeippen, C. J. 2000, MNRAS, 312, 813 

Tanaka, Y. & Shibazaki, N. 1996 ARA&A, 34, 607 

Tolea, A., Krolik, J. H. & Tsvetanov, Z. 2002, ApJ, 578, L31 

Turnshek, D. A., Monier, E. M., Sirola, C. J., & Espey, B. R. 1997, ApJ, 476, 40 

Veilleux, S. & Osterbrock, D. E. 1987, ApJS, 63, 295 

Vestergaard, M. 2002, ApJ, 571, 733 



Weymann, R.J., Morris, S.L., Foltz, C.B., k Hewett, P.C. 1991, ApJ, 373, 23 
Woo, J.-H. & Urry, C. M. 2002, ApJ, 579, 530 



This preprint was prepared with the AAS IATgX macros v5.0. 



